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ABSTRACT 

We  consider  several  problems  involving  points  and  planes  in 
three  dimensions.  Our  results  include:  (i)  The  maximum  number  of 
faces  bounding  m  distinct  cells  in  an  arrangement  of  n  planes  is 
O^m^'^n  log  n  +  n^)\  one  can  calculate  m  such  cells,  specified  by  a 
point  in  each,  in  worst-case  time  0{m^'^n  log-'n  +  n^log  n).  (ii)  The 
maximum  number  of  incidences  between  m  points  and  n  planes  is 
O {m^'^ n  log  n  +  n^),  provided  these  points  are  vertices  of  the 
arrangement  of  these  planes,  but  this  number  is  only 
O^m^'^'^n'^'^'^^^  +  m  +  n  log  m),  for  any  8>0,  for  any  collection  of 
points  no  three  of  which  are  collinear.  (iii)  For  an  arbitrary  collection 
of  m  points,  one  can  calculate  the  number  of  incidences  between  them 
and  n  planes  by  a  randomized  algorithm  whose  expected  time  com- 
plexity is  0(  im'^"^~^n^'^*^^  +  m)  log^n  +  n  log  n  logm),  for  any 
8>0.  (iv)  Given  m  points  and  n  planes,  one  can  find  the  plane  lying 
immediately  below  each  point  in  randomized  expected  time 
0{{m^'^~^n'^"^"''^^  +  m  ]  \og^n  +  n  log  n  logm),  for  any  8>0.  (v) 
The  maximum  number  of  facets  (i.e.  (J  — l)-dimensional  faces) 
bounding  m  distinct  cells  in  an  arrangement  of  n  hyperplanes  in  d 
dimensions,  d>2>,  is  0{m^'^n'^'^  log  n  +  n'^~^).  This  is  also  an  upper 
bound  for  the  number  of  incidences  between  m  points  and  n  hyper- 
planes in  d  dimensions,  provided  the  points  are  vertices  of  the 
arrangement  of  these  hyperplanes.  The  combinatorial  bounds  in  (i), 
(v)  and  the  general  bound  in  (ii)  are  almost  tight  in  the  worst  case. 
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posium, Tokyo,  Japan  1988. 
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1.   Introduction 

In  this  paper  we  extend  the  techniques  developed  in  a  companion  paper  [EGS] 
to  a  variety  of  problems  involving  arrangements  of  planes  in  three  dimensions  or 
hyperplanes  in  higher  dimensions,  as  listed  in  the  abstract.  Problems  of  this  kind,  in 
particular  those  involving  incidences  between  points  and  lines,  or  hyperplanes,  or 
various  other  types  of  curves  or  surfaces,  have  been  considered  for  a  long  time  as 
key  problems  in  combinatorial  geometry  (see  [PE]  for  a  survey).  A  major  break- 
through was  achieved  by  Szemeredi  and  Trotter  [ST]  who  have  shown  that  the  max- 
imum number  of  incidences  between  m  points  and  n  lines  in  the  plane  is 
0(m^'^n^'^  +  m  +  n) .  Edelsbrunner  and  Welzl  [EW]  have  extended  the  lower 
bound  to  also  apply  to  the  complexity  of  m  distinct  cells  in  an  arrangement  of  n  lines 
in  the  plane  (see  also  [Ca]),  and  have  also  obtained  a  weaker  upper  bound  for  that 
complexity.  These  lower  bound  and  weak  upper  bound  have  later  been  extended  to 
three  and  higher  dimensions  by  Edelsbrunner  and  Haussler  [EH]  and  by 
Edelsbrunner  [E],  showing  that  the  maximum  complexity  of  m  distinct  cells  in  an 
arrangement  of  n  hyperplanes  in  d^3  dimensions  is  flim^'^n'^'^  +  n^~^)  and 
Oim^'^n'^'^  +  n''~^);  here  the  complexity  of  a  cell  is  the  number  of  facets  (i.e. 
(^  — l)-dimensional  faces)  bounding  it.  The  lower  bound  is  derived  by  a  fairly  easy 
"lifting"  of  the  two-dimensional  construction  of  [ST],  [EW]  to  higher  dimensions. 

Recent  progress  has  been  achieved  by  the  authors  in  two  companion  papers 
[CEGSW],  [EGS].  These  papers  introduce  a  new  technique,  whose  aim  is  to  decom- 
pose problems,  such  as  the  above,  into  several  subproblems  of  smaller  size;  such  a 
decomposition  is  achieved  by  drawing  a  random  sample  either  of  the  given  points  or 
of  the  given  lines  or  curves,  and  then  splitting  the  collection  of  points  and  curves 
into  subcollections  according  to  their  interaction  with  the  sample  objects.  This 
decomposition  leads  to  a  recurrence  formula  for  the  desired  quantity  which  is  then 
solved  to  yield  improved  upper  bounds.  This  approach  has  led  to  a  tight  bound 
&(m^'^n^'^  +  n)  on  the  complexity  of  m  distinct  cells  in  an  arrangement  of  n  lines 
in  the  plane,  and  to  several  other  improved  upper  bounds  for  the  maximum  number 
of  incidences  between  points  and  curves,  as  well  as  for  the  maximum  complexity  of 
many  cells  in  arrangements  of  curves,  where  the  curves  are  either  line  segments,  or 
circles,  or  pseudolines. 

These  techniques  can  also  be  adapted  to  yield  comparably  efficient  randomized 
algorithms  for  calculating  m  cells  in  a  given  arrangement  or  for  finding  all 
incidences  between  given  points  and  curves  (see  e.g.  [EGHSSSW]). 

In  this  paper  we  extend  these  recent  techniques  to  three  and  higher  dimensions, 
to  obtain  the  results  stated  in  the  abstract.  The  bounds  that  we  obtain  (for  cases  (i), 
(v),  and  the  general  bound  in  (ii))  are  almost  tight  in  the  worst  case,  apart  for  a  log 
factor,  as  follows  from  [EH],  [E].  Our  techniques  rely  heavily  on  the  bounds 
obtained  for  the  corresponding  two-dimensional  problems.  In  certain  cases  (notably 
in  cases  (i)  and  (v))  our  decomposition  strategy  is  considerably  simpler  than  that 
used  in  the  two-dimensional  cases,  mainly  because  the  bounds  we  want  to  establish 
are  larger;  in  the  other  cases,  we  use  a  more  complex  decomposition  strategy  which 
is  a  generalization  of  the  2-d  strategies. 

The  paper  is  organized  as  follows.  In  section  2  we  analyze  the  complexity  of 
many  cells  in  arrangements  of  planes  in  three  dimensions.    Section  3  presents  an 
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efficient  worst-case  algorithm  for  their  calculation.  In  section  4  we  consider  the 
"simpler"  problems  (ii)-(iv)  for  three-dimensional  arrangements  of  planes,  and  in 
section  5  we  generalize  the  analysis  of  section  2  to  higher  dimensions. 

2.   The  complexity  of  many  cells  in  an  arrangement  of  planes 

Let  n  =  {iTi,  .  .  .  jiTn}  be  a  collection  of  n  planes  in  three  dimensions,  and  let 

ci c„  be  m  distinct  cells  of  the  arrangement  A=A(U)  of  these  planes.  (We 

assume  familiarity  of  the  reader  with  the  notion  of  arrangements;  see  [E]  for  basic 
information  concerning  arrangements.)  The  complexity  of  a  cell  is  the  number  of 
(two-dimensional)  faces  bounding  the  cell  (by  Euler's  formula,  this  complexity  also 
bounds,  up  to  a  constant  factor,  the  number  of  vertices  and  edges  of  the  cell).  Our 
goal  is  to  obtain  sharp  bounds  on  the  maximum  complexity  K{m,n)  of  m  such  cells. 
A  lower  bound  K{m,n)  =  n{m^'^n  +  n^)  has  been  established  in  [EH]  (see  also 
[E]).  Here  we  prove  an  almost  matching  upper  bound  of  the  form 
Kim,n)  =  Oim'^'^n  log  n  +  n'^). 

It  is  convenient  to  reformulate  the  problem  as  the  following  one.  Let 
Pi.  •  •  •  ,Pm  be  m  given  points.  Our  task  is  to  analyze  the  complexity  of  the  cells  of 
A  (11)  that  contain  these  points,  where  any  such  cell,  even  if  it  contains  more  than 
one  point,  is  to  be  counted  just  once.  We  will  however  assume  that  no  point  lies  on 
any  plane,  and  that  initially  no  two  of  these  points  lie  in  the  same  cell.  (The  problem 
is  formulated  in  a  more  general  way,  because  our  analysis  will  break  the  problem 
recursively  into  subproblems  where  cells  may  contain  more  than  one  point.) 

The  analysis  uses  the  following  straightforward  divide  and  conquer  approach. 
We  first  partition  11  into  two  subsets  IIi,  112,  each  containing  roughly  nil  planes;  we 
will  refer  to  the  planes  in  these  subsets  as  the  red  planes  and  the  blue  planes  respec- 
tively. We  next  calculate  recursively  the  cells  containing  the  points  /?,•  in  each  of  the 
subarrangements  A(ni),  A  (Ha),  and  then  intersect  these  cells  to  obtain  the  desired 
cells  of  A  (11).  Note  that  if  a  point  p  lies  in  a  red  cell  R  and  in  a  blue  cell  B,  then  the 
final  cell  E  containing  p  is  simply  Rr\B,  and,  since  all  these  cells  are  convex  polyhe- 
dra,  the  complexity  of  E  is  bounded  by  the  sum  of  the  complexities  of  R  and  of  B. 
However,  R  (or  B)  might  contain  more  than  one  point  p,-,  in  which  case  we  would 
not  want  to  use  the  entire  complexity  of  R  in  bounding  the  complexity  of  each  of  the 
final  cells  containing  these  points,  as  this  duplication  might  result  in  an  unacceptably 
high  bound.  This  major  technical  difficulty  is  dealt  with  in  the  following  subsection, 
and  its  solution  (the  "combination  lemma")  then  enables  us  to  complete  our  analysis 
o^K{m,n). 

2.1.   The  combination  lemma  for  planes  in  3-d  space 

Let  pi,  .  .  .  ,p„  be  m  points.  We  are  given  two  arrangements  of  planes,  the 
red  arrangement,  consisting  of  rir  red  planes,  and  the  blue  arrangement,  consisting 
of  rib  blue  planes.  Let  /?i,  .  .  .  .^^  be  the  cells  in  the  red  arrangement  containing  the 
points  Pi  (where  a  cell  can  contain  more  than  one  point),  and  let  fli,  .  .  .  ,S,  be  the 
blue  cells  containing  these  points.  For  each  p,-  let  £,■  =  R^^  n  5,,  be  the  "purple" 
cell  containing  p,-  (where  R^.,,  B,.  are  the  red  and  blue  cells  containing  that  point). 
By  considering  only  a  subset  of  the  points  p,-,  if  necessary,  we  can  assume,  without 
loss  of  generality,  that  all  the  cells  £,•  are  distinct,  and  thus  have  pairwise  disjoint 
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interiors.  (While  this  assumption  explicitly  holds  at  the  first  level  of  recursion, 
further  levels  will  generally  have  more  than  one  point  in  the  cells  £,;  to  enforce  the 
above  property,  we  simply  choose  one  representative  point  out  of  each  cell  E,  and 
(temporarily)  ignore  the  others.)  Let  p,  P  denote  the  total  number  of  faces  on  the 
boundaries  of  the  red  cells  and  the  blue  cells  respectively.  Our  goal  is  to  obtain  a 
sharp  upper  bound,  in  terms  of  p  and  P,  on  the  total  number  of  faces  on  the  boun- 
daries of  the  purple  cells  E,-. 

As  in  the  2-d  case  [EGS]  and  as  noted  above,  it  is  easily  verified  that  if  each 
red  cell  and  each  blue  cell  contains  just  a  single  point,  then  the  desired  purple  com- 
plexity is  at  most  p  +  p.  The  main  step  in  bounding  the  purple  complexity  will  thus 
aim  to  partition  each  red  cell  that  contains  more  than  one  point  into  subcells,  each 
containing  only  a  single  point,  using  (portions  of)  the  blue  planes,  and  to  obtain  a 
symmetric  partitioning  of  the  blue  cells. 

We  consider  the  blue  planes  in  a  fixed  sequence,  and  add  them  to  the  red 
arrangement  one  at  a  time.  We  use  each  of  them  to  split  as  many  red  cells  as  neces- 
sary into  subcells,  called  "red-blue"  cells,  as  follows.  Let  R  =  /?,  be  a  fixed  red  cell 
which  contains  m,>l  points.  Take  the  first  blue  plane  ir  in  the  sequence  which 
passes  through  R  and  splits  the  points  in  R  into  two  non-empty  subsets  (our  assump- 
tions imply  that  there  exists  at  least  one  such  plane).  The  plane  it  partitions  R  into 
two  subcells.  If  at  least  one  of  them  still  contains  more  than  one  point,  we  take  a 
subsequent  blue  plane  which  splits  the  points  in  that  subcell  into  two  nonempty  sub- 
sets, and  partition  that  subcell  further  into  two  subcells  by  this  plane.  We  continue 
this  process  until  each  resulting  red-blue  subcell  contains  just  one  point.  This  pro- 
cess is  applied  simultaneously  to  all  other  red  cells. 

Some  features  of  this  splitting  process  are  worth  noting.  A  blue  plane  it  can 
(and  will)  be  used  to  split  more  than  one  subcell.  That  is,  it  will  split  into  two  sub- 
cells  each  red-blue  cell  that  it  crosses,  provided  it  also  splits  the  set  of  points  within 
that  subcell  into  two  nonempty  subsets.  However,  if  ir  crosses  a  red-blue  cell 
without  splitting  the  set  of  points  within  it  (e.g.  if  that  cell  contains  just  a  single 
point),  then  the  cell  remains  unchanged. 

What  we  want  to  estimate  is  the  "red  complexity"  of  the  red-blue  cells,  namely 
the  number  of  red  (sub)faces  bounding  these  cells  (so  the  blue  faces  along  their 
boundaries  are  ignored);  more  precisely,  we  want  to  estimate  the  number  of  addi- 
tional red  subfaces  generated  by  the  intersections  of  each  blue  plane  with  the 
current  set  of  red-blue  cells. 

Consider  the  cell  splitting  caused  by  the  i-th  blue  plane  it,.  Let  Rj  be  a  red  cell 
that  has  already  been  split  into  a  number  of  red-blue  subcells,  and  suppose  that  it,- 
splits  k  =  kj  of  these  subcells  (in  the  sense  of  also  splitting  the  sets  of  points  they 
contain).  Since  the  red-blue  subcells  are  convex  polyhedra  by  construction,  the 
intersection  of  it,-  with  each  of  them  is  a  convex  polygon  q  (where  these  polygons 
have  pairwise  disjoint  interiors),  contained  in  the  intersection  c  =  Cj  =  'ni  D  Rj  which 
itself  is  a  convex  polygon.  To  each  red  face  /of  a  red-blue  subcell  which  is  split  into 
two  subfaces  by  it,-,  there  corresponds  a  subsegment  of  some  edge  of  c  (which  is  the 
intersection  of /with  it,).  Thus  it  follows  easily  that  the  total  increase  Uj  in  the  red 
complexity  of  the  red-blue  subcells  within  Rj  as  caused  by  it,-  is  equal  to  the  number 
of  "red"  edges  of  the  red-blue  intersection  polygons  q,  namely  edges  that  he  along 
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the  boundary  of  c  (see  Figure  1).  But,  as  in  the  proof  of  the  Combination  Lemma  1 
of  [EGS],  the  total  number  of  red  edges  of  k  convex  polygons  (with  disjoint  interi- 
ors) contained  in  a  red  convex  polygon  c,  is  at  most  \c  \+2k  —  2,  where  \c  \  denotes 
the  number  of  edges  of  c  (and  where  a  red  edge  means  an  edge  that  lies  on  the 
boundary  of  c). 

Hence,  summing  these  increases  /y  over  all  red  cells  Rj,  it  follows  that  the  total 
increase  in  the  red  complexity  caused  by  ir,-  is  at  most  ^  (  \cj\  +  2kj  —  2),  where  the 

;' 
summation  extends  over  all  red  cells  Rj  containing  at  least  one  red-blue  subcell  that 

has  been  split  by  it,-;  we  let  ^j  denote  the  number  of  such  red  cells. 

Consider  first  the  sum  ^  \cj  \.  Let  A  be  the  arrangement  of  the  rir  lines  formed 

by  intersecting  it,-  with  each  of  the  rir  red  planes.  Then  our  sum  gives  the  overall 
complexity  of  ^,  faces  (i.e.  the  polygons  Cj)  of  A.  By  the  results  of  [CEGSW],  this 
complexity  is  at  most  O^^j'^n}'^  +  rir).  Summing  this  over  all  blue  planes  it,-,  we 
obtain 

2  0(?P/i2/3  +  rir)  =  0((2  ^V^)-nV^  +  nrrib)  . 
.=1  1=1 

But 

rtb  rib 

V    c2/3  <:  /  V    C  ^2/3.„l/3  <  „2/3_l/3 

«=1  1=1 

Because  2  ^i  gives  the  number  of  times  a  red  cell  has  been  split  by  a  blue  plane 
(namely  a  red-blue  subcell  within  it  has  been  so  split),  and  this  number  is  clearly 
bounded  by  the  total  number  of  points  m. 

Finally  note  that  2  ^j  's  the  total  number  of  subcell  splittings  caused  by  it,-. 

Since,  as  just  noted,  the  overall  number  of  subcell  splittings  is  at  most  m,  it  follows 
that  the  sum  of  the  corresponding  terms  2  (2kj-2)  over  all  blue  planes  is  at  most 

2m.  Hence  the  total  increase  in  the  red  complexity  is  at  most 

Oim^'^nl'^n^'^  +  rirTib  +  m)  . 

We  now  obtain  in  a  completely  symmetric  manner  a  partitioning  of  the  blue 
cells  into  a  collection  of  "blue-red"  subcells,  each  containing  just  a  single  point,  and 
derive  a  symmetric  bound  on  their  overall  "blue"  complexity.  Note  that  if  a  point  p,- 
lies  in  a  red-blue  subcell  c,-  and  in  a  blue-red  subcell  c,-,  then  the  final  purple  cell  f,- 
containing  p,-  is  E,-  =  c,-  PI  c,.  This  is  because  £",-  =  /J^,.  fl  B,.  by  definition,  and  by 
construction  £,•  is  contained  in  both  c,-  C  /J^,.,  c,-  C  fl,,.  Moreover,  each  red  face  of 
Ei  must  be  contained  in  a  red  face  of  c,-,  and  each  blue  face  of  E,  must  be  contained 
in  a  blue  face  of  c,-.  Consequently,  the  number  of  faces  of  £",  is  bounded  by  the  red 
complexity  of  c,-  plus  the  blue  complexity  of  c,-. 

Putting  everything  together,  we  finally  obtain 

Lemma  2.1  (Combination  Lemma  for  Arrangements  of  Planes):  Given  an  arrange- 
ment of  rir  red  planes,  another  arrangement  of  rih  blue  planes,  and  m  points  not 
lying  on  any  of  these  planes,  the  total  number  of  faces  bounding  the  "purple"  cells 


-6 


containing  these  points,  in  the  arrangement  formed  by  all  rir  +  ni,  planes,  is  at  most 

p   +   P  +  0[m2'3(ni'3„2/3  +  „2/3„l/3)   +  „^„^  +  ^]   ^ 

where  p  (resp.  p)  is  the  total  number  of  faces  bounding  the  cells  of  the  red  (resp. 
blue)  arrangement  containing  these  points. 


2.2.   The  Complexity  of  Many  Cells 

We  now  return  to  our  analysis  of  K(m,n).  Applying  the  Combination  Lemma 
to  the  two  subarrangements  A(IIi),  AClli),  we  can  readily  obtain  a  recurrence  of 
the  form 

K(m,n)  =  2K(m,n/2)  +  0(m^'^n  +  n^  +  m)  . 

However,  we  do  not  continue  this  recursive  process  all  the  way  down.  When  m 
becomes  greater  than  or  equal  to  n^  we  use  the  trivial  bound 
K(m,n)  =  0{n^)  =  0{m).  Since  when  this  happens  we  have  m^(2n)3,  it  follows 
that  m  =  m^'^m^'^  =  0{m^'^n).  By  the  same  argument,  we  can  drop  the  term 
0(m)  in  the  above  recurrence,  since  we  apply  it  only  for  m^n  ,  in  which  case  it  is 
dominated  by  O(m^'^n).  It  is  now  easy  to  verify  that  the  solution  to  the  above 
recurrence  is 

K(m,n)  =  Oim'^'^n  log  n  +  n^)  . 

(Note  that  the  log  enters  the  formula  for  K{m,n)  because  in  the  recurrence  relation 
the  term  O(m^'^n)  is  linear  in  n.) 

We  have  thus  shown 

Theorem  2.2:  The  complexity  of  m  distinct  cells  in  an  arrangement  of  n  planes  in 
three  dimensions  is  0{m^'^n  log  n  +  n^). 

Corollary  2.3:  The  maximum  number  of  incidences  between  m  points  and  n  planes 
in  three  dimensions,  where  it  is  assumed  that  each  point  is  a  vertex  of  the  arrange- 
ment of  these  planes,  is  Oim^'^n  log  n  +  n^). 

Proof:  We  use  the  following  well-known  argument  (see  [EW],  [E]).  Replace  each 
plane  by  two  parallel  planes  translated  a  small  distance  away  in  both  directions.  If 
all  these  translation  distances  are  equal  and  sufficiently  small,  each  vertex  p  of  the 
arrangement  of  the  original  planes  will  lie  in  a  distinct  cell  of  the  new  arrangement, 
whose  complexity  is  2d,  where  d  is  the  number  of  original  planes  incident  to  p.  Thus 
the  maximum  number  of  incidences  between  m  points  and  n  planes  as  above  is 

bounded  by  —Kim,  2n),  which  completes  the  proof  of  the  corollary.  □ 

Remarks: 

(1)  The  proof  of  Theorem  2.2  is  significantly  simpler  than  the  elaborate  proof  of  the 
corresponding  bound  in  two  dimensions  (see  [CEGSW],  [EGS]),  and  of  the  simpler 
3-d  variants  analyzed  in  the  following  section.  This  is  because  (i)  we  have  made 
strong  use  of  the  2-d  results  in  the  proof  of  the  combination  lemma,  and  (ii)  the 
weaker  bound  in  three  dimensions  allows  a  straightforward  divide  and  conquer  stra- 
tegy, which  cannot  be  used  to  derive  the  sharper  2-d  bound  (or  the  sharper  3-d 


bounds  for  the  problems  given  below). 

(2)  The  bound  given  in  Corollary  2.3  almost  matches  the  lower  bound 
n(m^'-'n  +  n^)  given  in  [EH],  [E].  It  is  interesting  to  note  that  the  lower  bound 
construction  in  [EH]  is  a  fairly  trivial  extension  of  the  2-d  lower  bound  for 
incidences  between  points  and  lines.  Thus  the  main  tools  used  in  the  derivations  of 
both  the  lower  and  upper  bounds  for  the  case  of  planes  are  principally  based  on 
results  on  2-d  arrangements. 

(3)  One  could  also  give  a  direct  proof  of  Corollary  2.3,  without  reducing  incidences 
to  the  complexity  of  many  cells.  One  such  direct  proof  is  an  easy  modification  of  the 
proof  given  above,  including  a  modified  combination  lemma  for  "merging 
incidences"  between  a  set  of  points  and  two  collections  of  red  and  of  blue  planes. 
As  a  hint  to  the  interested  reader,  who  may  wish  to  work-out  the  details  of  this 
modified  proof,  we  note  that  the  main  technical  difficulty  that  arises  is  that  the  given 
points  are  assumed  to  be  vertices  of  the  overall  arrangement,  but  they  need  not  be 
vertices  of  the  blue  or  of  the  red  subarrangements.  A  combination  lemma  is  then 
needed  to  handle  points  of  this  sort  (which  generally  lie  on  an  edge  of  one  arrange- 
ment and  on  a  plane  of  the  other). 

(4)  An  obvious  open  problem  is  whether  the  upper  bounds  in  Theorem  2.2  or  in 
Corollary  2.3  can  be  tightened  to  match  the  corresponding  lower  bounds,  i.e. 
whether  the  log  n  factor  can  be  removed.  This  factor  is  just  an  artifact  of  our 
simple-minded  divide-and-conquer  strategy,  and  one  would  hope  that  better  stra- 
tegies might  get  rid  of  this  factor. 

(5)  The  technique  presented  in  this  section  has  recently  been  generalized  (and  some- 
what simplified)  in  [AS]  to  obtain  sharp  upper  bounds  on  the  complexity  of  cells  in 
an  arrangement  of  n  (possibly  intersecting)  triangles  in  3-space.  For  example,  it  is 
shown  there  that  the  total  complexity  of  all  non-convex  cells  in  such  an  arrangement 
is  0(/i^'-'a(«)^'-'log'*'^/i),  which  is  abnost  optimal  in  the  worst  case. 

(6)  It  would  be  interesting  to  extend  these  results  to  arrangements  of  other  surfaces 
(e.g.  spheres)  in  three  dimensions.  Since  the  analysis  given  above  strongly  uses  the 
straightness  of  planes  and  convexity,  new  techniques  will  be  required  to  handle 
curved  surfaces. 

3.   Calculating  Many  Cells 

Next  we  sketch  an  algorithm  for  efficient  calculation  of  m  cells  in  an  arrange- 
ment of  n  planes  in  time  0(m^'^n  log^n  +  n^log  n)  in  the  worst  case. 

The  algorithm  relies  on  the  following  operation.  Suppose  we  are  given  two 
intersecting  convex  polyhedra  P  and  Q  in  three  dimensions,  of  size  p  and  q  respec- 
tively. We  are  told  that  P  and  Q  intersect  and  we  are  furthermore  given  a  point  x  in 
their  intersection.  Our  goal  is  to  construct  R  =  Pr\Q  in  time  proportional  to  r,  the 
size  of  ^,  times  a  polylogarithmic  factor  in  p,  q,  and  r.  (We  think  of  p  as  represent- 
ing the  number  of  faces  of  P,  which  is  linearly  related  to  the  number  of  edges  and 
the  number  of  vertices  of  P;  similarly  for  q  and  r.) 

In  order  to  do  this  we  will  first  develop  a  subroutine  that  intersects  a  planar 
convex  polygon  with  a  convex  polyhedron  in  three  dimensions  within  an  analogous 
time  bound.  So  let  C  be  a  convex  polygon  of  size  c,  and  let  Q  and  q  be  as  above. 


Again,  we  assume  that  we  are  given  a  point  z  in  the  intersection  5  of  C  and  Q, 
where  5  is  a  planar  convex  polygon  of  size  s. 

It  is  actually  instructive  to  go  down  in  dimension  one  more  time  and  consider 
how  to  intersect  two  convex  polygons  in  output-sensitive  time.  A  straightforward 
technique  for  doing  so  is  described  in  [EGS]  and  is  based  on  "shooting".  The  idea 
is  that  we  find  a  point  on  the  boundary  of  the  intersection  and  then  walk  along  this 
boundary  by  shooting  from  the  current  edge  of  the  inner  polygon  towards  the  outer 
one,  in  order  to  find  the  next  vertex  of  the  intersection.  There  are  two  cases.  It  may 
happen  that  the  inner  polygon  edge  finishes  before  it  exits  the  outer  polygon,  in 
which  case  we  move  to  the  next  edge  of  the  inner  polygon  and  repeat  the  process. 
Or  we  may  encounter  the  outer  polygon  before  the  inner  polygon  edge  finishes,  in 
which  case  we  have  a  ne\y  vertex  of  the  intersection  polygon  and  the  roles  of  the 
inner  and  outer  polygons  are  swapped.  Since  shooting  a  ray  in  a  convex  polygon 
takes  logarithmic  time,  the  whole  intersection  procedure  can  be  done  in  time  pro- 
portional to  the  size  of  the  intersection  times  a  logarithmic  factor.  Note  that  loga- 
rithmic shooting  within  a  convex  polygon  can  be  done  without  any  preprocessing. 

We  now  discuss  how  to  "lift"  this  argument  to  three  dimensions  and  solve  our 
subproblem.  Note  first  of  all  that  ray  shooting  in  a  convex  polyhedron  can  also  be 
done  in  logarithmic  time,  for  example  by  using  a  Dobkin-Kirkpatrick  style  [DK] 
hierarchical  representation  (see  also  [E,  Section  9.5]).  So  we  proceed  as  follows. 
From  our  given  point  z  in  the  intersection  of  C  and  Q  shoot  any  ray  in  the  plane  of 
C  and  find  its  nearest  intersection  point  v  with  C  or  Q.  We  now  implement  a 
"walk-around"  procedure  similar  to  the  one  outlined  above.  If  v  is  on  C,  then  we 
shoot  along  the  current  edge  of  C  to  see  if  we  exit  Q.  As  long  as  have  not,  we  just 
keep  walking  around  C.  If  we  do  exit  Q  then  we  have  a  new  vertex  w  of  5,  lying  in 
some  face  f  oi  Q.  Now  things  are  more  interesting,  because  we  need  to  shoot  from 
w  in  both  C  and  the  face  f  of  Q  in  order  to  find  the  next  vertex  of  S.  It  is  still  true, 
however,  that  both  of  these  shots  can  be  accomplished  in  logarithmic  time,  and  this 
remains  so  as  we  go  around  along  the  surface  of  Q.  So  the  polygon-polyhedron  case 
is  only  slightly  more  complex  than  the  polygon-polygon  case  .  The  intersection  S  can 
again  be  computed  in  time  proportional  to  s  times  a  logarithmic  factor. 

We  now  go  back  to  the  original  polyhedron-polyhedron  problem.  From  the 
common  point  x  oi  P  and  Q  let  us  shoot  an  arbitrary  ray  and  let  the  nearest  face  hit 
be  C,  say  it  is  a  face  of  P.  Then  we  can  apply  the  above  procedure  to  intersect  C 
with  Q  and  obtain  S,  a  face  of  the  desired  intersection  output  R.  Now  let  e  be  an 
edge  of  S.  We  claim  that  our  subprocedure  above  can  be  used  to  compute  the  face  T 
of  R  adjacent  to  5  along  e.  There  are  three  cases.  If  g  is  a  subedge  of  an  edge  of  P, 
then  we  simply  take  the  other  face  of  P  along  that  edge  and  intersect  it  with  Q  to  get 
T.  If  g  is  a  subedge  of  an  edge  of  Q  we  do  the  symmetric  operation.  Finally  if  e  is  an 
edge  arising  out  of  the  intersection  of  a  face  (in  this  case  C)  of  P  with  a  face  D  of 
Q,  then  we  simply  switch  the  roles  of  P  and  Q  (inner  with  outer)  and  now  intersect 
polygon  D  with  polyhedron  P.  If  we  repeat  this  process  we  will  eventually  discover 
all  the  faces  of  the  intersection  polyhedron  R.  It  is  clear  that  the  cost  of  the  compu- 
tation is  proportional  to  the  number  of  edges  of  R,  times  a  factor  which  is  the  cost 
of  shooting,  paid  a  constant  number  of  times  per  edge.  Overall  this  gives  us  the 
desired  output  sensitive  algorithm  for  polyhedral  intersection. 
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We  now  apply  this  technique  to  compute  the  m  desired  cells  in  an  n  plane 
arrangement.  The  cells  are  assumed  to  be  specified  by  m  points  pi,  .  .  .  ,pm,  one 
point  lying  in  the  interior  of  each  cell.  We  apply  a  divide  and  conquer  technique  as 
in  the  preceding  proof.  That  is,  we  divide  the  collection  of  planes  into  two  subcoUec- 
tions  ("blue"  and  "red"),  each  containing  roughly  n/2  planes,  and  recursively  com- 
pute the  cells  containing  the  given  points  in  the  blue  and  in  the  red  subarrange- 
ments.  We  next  preprocess  each  of  these  convex  polyhedral  cells,  as  in  [DK],  to 
prepare  for  logarithmic-time  ray  shooting  inside  them.  The  total  cost  of  preprocess- 
ing is  bounded  by  0{K(m,nl2)).  Then,  for  each  point  pj,  we  apply  the  above 
polyhedron-polyhedron  intersection  routine.  However,  since  two  of  these  points 
could  lie  in  the  same  blue  and  the  same  red  cell  (at  least  when  we  are  deep  down 
the  recursive  process),  we  have  to  exercise  some  care  to  make  sure  that  no  intersec- 
tion polyhedron  is  computed  more  than  once.  This  however  is  easy  to  enforce  by 
checking  for  each  newly  computed  vertex  of  an  intersection  polyhedron  that  it  has 
not  been  previously  obtained  as  a  vertex  of  another  intersection.  With  this  precau- 
tion, the  output-sensitive  behavior  of  the  above  routine  implies  that  all  desired  cells 
can  be  computed  in  time  0(K(m,n)  log  n)  (including  the  preprocessing  overhead  as 
well).  Moreover,  when  m^n^ ,  we  can  simply  calculate  the  entire  arrangement  in 
time  0(n^)  =  0(m)  (as  in  [EOS]),  and  pick  out  the  desired  cells.  To  accomplish  this 
last  step,  we  apply  the  point  location  algorithm  in  [EOS]  (or  that  in  [Ch])  to  locate 
each  of  the  given  m  points  in  the  computed  arrangement,  in  total  time  0(m  log^n). 
All  these  considerations  lead  to  the  following  recurrence 

Tim  n)  =  |27'('«.«/2)  +  0(m^'\  log^n  +  n^\og  n)  ,     m<n^ 
{0(m  log^rt)  ,  m>n^ 

for  the  maximum  time  T(m,n)  required  by  the  algorithm;  this  recurrence  is  easily 
seen  to  solve  to 

T{m,n)  =  0(m^'^n  log^n  +  n^log  n  +  m  log^/i)  . 

Finally,  if  we  wish  to  calculate  m  distinct  cells  in  an  arrangement  of  planes,  and  the 
cells  are  given  by  a  single  point  in  each,  then  necessarily  m<n^ ,  so  that  we  can  drop 
the  term  0(m  log^n)  because  it  is  then  dominated  by  the  first  term  of  this  bound. 
However,  when  we  are  not  guaranteed  that  each  desired  cell  contains  a  unique 
designating  point,  the  last  term  may  become  dominant.  We  thus  conclude 

Theorem  2.4:  One  can  calculate  m  distinct  cells  in  an  arrangement  of  n  planes  in  3- 
space,  when  the  cells  are  designated  by  a  unique  point  in  each,  in  time 
0{m^'^n  log^/i  +  /i^log  n)  in  the  worst  case. 

4.   Simpler  problems  involving  points  and  planes 

Let  n  =  {iTi,  .  .  .  ,'n'n}  be  a  collection  of  n  planes  in  3-dimensional  space,  and 
let  pi,  .  .  .  ,p„  be  m  points.  In  this  section  we  consider  several  variants  of  the 
many-cells  problem  studied  above,  involving  these  points  and  planes.  Specifically, 
we  consider  the  following  problems: 

(a)  Determine  for  each  point  p,-  the  plane  itj  lying  immediately  below  it. 

(b)  Count  the  number  of  incidences  between  the  given  points  and  planes. 
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(c)    Obtain  an  upper  bound  on  the  number  of  incidences  between  these  points  and 

planes,  assuming  that  no  three  points  are  coUinear  (or,  symmetrically,  that  no 

three  planes  pass  through  a  common  line). 
Note  that  problem  (a),  and,  with  some  extra  care,  problem  (b)  as  well,  can  be 
solved  by  the  algorithm  of  section  2.3.  However,  as  shown  below,  these  variants  (as 
well  as  problem  (c))  admit  solutions  with  bounds  that  are  significantly  lower  than 
those  obtained  in  the  previous  section. 

We  will  present  in  detail  a  solution  to  the  first  problem,  comment  briefly  on  the 
modifications  necessary  to  obtain  a  similar  solution  to  the  second  problem,  and  then 
analyze  the  third  problem,  which  happens  to  have  a  better  bound  than  the  others. 

We  will  use  an  approach  similar  to  that  used  in  [EGS]  for  problems  involving 
two-dimensional  arrangements;  this  approach  involves  a  divide-and-conquer  strategy 
based  on  random  sampling.  Specifically,  we  pass  to  the  dual  space  (using  standard 
duality  that  preserves  above/below  relationships;  see  [E]),  so  that  each  plane  ir,-  is 
transformed  to  a  point  irf ,  and  each  point  p^  becomes  a  "query"  plane  pf .  We  con- 
struct a  "customized"  partition  tree  T  for  the  set  of  n  dual  points  irf  and  for  the  col- 
lection of  m  dual  query  planes  pj .  Each  node  v  of  7  is  associated  with  a  subset  T\.* 
of  Tiv  of  these  dual  points,  and  with  a  collection  P?  of  m^  dual  planes  that  "reach"  v. 
At  each  node  v  for  which  rriy  "^  nl,  we  stop  the  recursive  construction  of  T.  Other- 
wise we  choose  a  random  sample  of  r  dual  planes  pf  (for  some  large  but  fixed 
integer  r),  form  their  arrangement,  and  decompose  its  cells  into  tetrahedra  with 
pairwise  disjoint  interiors;  the  number  of  tetrahedra  is  0{r^).  For  each  resulting 
tetrahedron  t  we  construct  a  child  w  =  Wt  of  v  in  T,  and  associate  with  w  the  set  II* 
of  all  dual  points  lying  in  t  and  the  set  P*  of  all  dual  planes  cutting  t  (these  are  the 
points  and  planes  that  are  said  to  "reach"  w).  (If  some  dual  point  irf  lies  on  a 
tetrahedron  boundary,  we  can  assign  it  to  any  of  the  adjacent  tetrahedra  without 
affecting  the  behavior  of  the  algorithm.)    The  e-net  theory  of  [HW]  or  the  similar 

probabilistic  analysis  of  Clarkson  [CI],  implies  that,  with  high  probability,  every 

cm 
such  tetrahedron  will  be  cut  by  at  most log  r  dual  planes,  for  some  positive  con- 
stant c. 

For  problem  (a)  we  store,  for  each  node  v  and  for  each  point  p,  whose  dual 
plane  reaches  v,  the  current  candidate  for  the  plane  ttj  lying  immediately  below  p,-. 
The  recursive  algorithm  can  then  proceed  as  follows. 

(i)  At  a  node  v  lying  at  the  bottom  of  the  recursion,  we  construct  the  full  arrange- 
ment of  the  set  Ilv  of  the  n^  planes  t7,-  associated  with  v,  and  locate  in  it  each  of  the 
wiv  corresponding  points  pj,  using  the  point  location  algorithm  of  [EOS]  or  [Ch]. 
This  can  be  accomplished  in  time  Oim^  log^/iv),  and  yields  initial  candidates  for  the 
planes  lying  below  each  of  these  points. 

(ii)  At  each  internal  node  v  we  process  the  set  B^  of  all  points  whose  dual  planes 
have  reached  the  parent  of  v  but  have  missed  v,  as  follows.  Let  p,  be  a  point  in  B^. 
Then  the  plane  pf  passes  either  above  all  the  points  ttj  €  Wt  or  below  all  of  them. 
By  duality,  the  point  p.-  lies  either  above  all  the  planes  tTj  ^  Ilv  or  below  all  of  them. 
In  the  latter  case,  no  plane  in  Ilv  can  be  the  plane  lying  directly  below  p,-,  so  these 
points  need  not  be  further  processed  at  v.  Points  p,-  of  the  first  type  lie  in  the  upper 
unbounded  cell  c  of  the  arrangement  4(11^.),  and  the  plane  of  Ilv   lying  directly 
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below 

Pi  is  the  one  containing  the  face  of  c  lying  directly  below  p,-.  We  thus  calculate  the 
cell  c  in  0(nv  log  n^)  time  (as  in  [PM]),  and  then  project  the  boundary  of  c,  as  well 
as  all  points  p,  €  B^  of  the  first  type,  onto  the  j^'-plane.  We  then  locate  each  of 
these  projected  points  in  the  planar  map  obtained  by  vertically  projecting  the  boun- 
dary of  c,  thereby  obtaining  the  plane  of  Ilv  lying  directly  below  p,-.  Thus  this  step 
can  be  accompHshed  in  (9((ibv  +  «v)  log  riy)  time,  where  by  is  the  size  of  fiy 
(iii)  Finally,  the  merge  step  at  each  internal  node  v  of  the  tree  is  performed  by  tak- 
ing, for  each  point  p,  reaching  v,  the  highest  candidate  plane  among  those  obtained 
at  each  of  the  children  of  v,  where  for  each  child  Wj  this  plane  is  either  obtained 
from  the  recursive  processing  of  wj  (if  p,  reaches  Wj),  or  from  processing  the  set 
Bwj  (if  Pi  belongs  to  that  set)  as  in  (ii)  above.  This  merge,  including  the  processing 
of  the  sets  fln-  for  all  children  w  of  v  (as  in  (ii)  above),  can  clearly  be  performed  in 
total  time  0((mv  +  nv)  log  «v)- 

Let  T(m,n)  denote  the  maximum  time  required  by  the  algorithm  for  an  input  of 
m  points  and  n  planes,  under  the  assumption  that  each  random  sampling  step  (at 
each  node  of  T)  does  indeed  produce  a  good  partition  of  the  dual  space  (see  below 
for  a  discussion  of  this  issue).   Taking  into  account  the  overhead  needed  to  construct 

thus  obtain  the  following  recurrence  formula  for  (T(m,n). 

'amlog'n  if  m  ^  n' 

M 

2  T{mi,ni)  +  {bm  +  b'n)  log  n      if  m  <  n^ 
j  =  l 


the  partition  tree  T,  we  t 
T{m,n)  < 


for  some  constants  a,  b,  b'  >  0,  where,  by  the  preceding  arguments,  the  m,,  «,  and 
M  satisfy  the  following  three  conditions 

(a)    M  =  0(r3); 

(b)   2  «.•  ^  «; 

1  =  1 
(c)    m,:£ log  r  for  each  i,  for  some  constant  c>0. 

Before  stating  our  results,  we  need  to  discuss  the  probabilistic  aspects  of  the 
above  algorithm.  Even  though  there  is  high  probability  of  choosing  a  good  random 
sample  at  any  single  node  v  of  T,  the  accumulated  probability  of  successful  drawings 
at  all  nodes  of  T  can  become  small,  unless  the  sample  size  increases  from  level  to 
level,  eventually  becoming  non-constant  near  the  bottom  of  recursion  and  thereby 
increasing  the  overhead  of  processing  a  single  node.  An  alternative  approach,  which 
is  the  one  we  follow,  is  to  verify,  at  each  step  of  the  algorithm,  that  the  currently 
sampled  subset  of  dual  planes  pf  is  indeed  a  good  sample,  namely  that  no  cell  in  its 

triangulated  arrangement  is  cut  by  more  than  0{—  log  r)  dual  planes.  This  can  be 

easily  verified  in  linear  time,  assuming  r  to  be  a  constant.  Thus,  with  no  increase  in 
the  asymptotic  complexity,  we  can  verify  that  the  sample  is  good.  If  the  sample  is 
bad,  we  simply  discard  it  and  choose  another.  Again,  the  results  of  [CI],  [HW] 
imply  that  the  expected  number  of  draws  until  a  good  sample  is  obtained  is  constant. 
Employing  this  approach,  we  obtain  a  randomized  algorithm,  whose  expected  time 
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complexity  is  proportional  to  T{m,n)  defined  above,  and  which  has  the  property 
that  after  a  node  is  processed,  the  resulting  sample  and  corresponding  space  parti- 
tioning are  good. 

Returning  to  the  analysis,  we  have 

Theorem  4.1:  One  can  determine  which  of  n  given  planes  lies  directly  below  each  of 
m  given  points,  using  a  randomized  algorithm  whose  expected  running  time  is 
bounded  by  [  Dm^^^'^n^"*"^^^  +  Am  ]  log^n  +  Bn  log  n  log  m,  for  any  8>0,  where 
the  coefficients  A,  B,  D  depend  on  S.  (Here  the  expectation  is  over  the  randomiza- 
tions performed  by  the  algorithm,  and  does  not  depend  on  the  input.) 
Proof:  If  we  follow  the  recurrence  given  above  for  T{m,n)  through  the  partition 
tree  T,  the  contribution  of  the  rightmost  term  b'riy  log  ny,  over  all  nodes  v  at  the 
same  recursive  level,  is  at  most  b'n  log  n  (see  the  analysis  in  [EGS]  for  a  similar 
argument).  Since  the  tree  T  has  only  log  m  levels  of  recursion,  it  follows  that  the 
overall  contribution  of  these  terms  is  0(/i  log  n  log  m).  It  is  therefore  sufficient  to 
drop  this  term  from  the  recurrence  defining  T{m,n),  and  prove  that  the  solution  T 
of  the  modified  recurrence  satisfies  T{m,n)  <  [  Dm^'^'^n^'^'^^^  +  Am  ]  log^n,  for 
any  8>0. 

Fix  8  >  0  and  choose  r  =  r(8)   to  be   sufficiently   large   (how  large   will  be 
apparent  from  subsequent  analysis). 

If  m  ^  n^  then  T(m,n)  ^  am  log^n  plainly  satisfies  the  required  inequality, 
assuming  A  ^  a.  So  suppose  m  <  n^ .  In  this  case 

^    _  •„ 3/4-6. _,  1/4+8   <;  „3/4-8_3/4  +  38  /•♦\ 

m  —  m  m  s  m  n  .  \  } 

By  induction  hypothesis  we  then  have 

T{m,n)  <  2  [Omj^^-^nV"^^^^  +  Ami  ]  log^/r,-  +  bm  log  n  . 
1  =  1 


M 

But  2  rrti  log^/i,  < 


Mc  ,        ^ 
log  r 

r 


■m  log^rt,  thus 


1 

,3/4-8„3/4+38 


^       -   •    -     -   ■    -  Mc 

r 


T(m,n)  ^  [  D  2  mi"*-''nf"'^'''  +  {A^  log  r  +  b)m  ]  log^n 
1 


Mc 
Thus,  using  (*)  and  putting  d  =  A log  r  +  b,  wc  obtain 

nm,n)  ^  [D  f  m^'^-^nf^^^s  +  ^^3/4-8„3/4+38  ^  ^^^2„ 

1 

-  t  ^^1^  \og^"-^r  m3"*-8.f  „?'4*38    +  ^^3/4-8^3/4+38  ^  ^^^2^ 
^  1 


But  since 

f  n^'^^^S   <  (^  „.)3/4  +  38.^1/4-38   <  y^^l/4-38  .^3/4  +  38 
1  1 

and  M  ^  c'r^  for  some  constant  c' ,  we  obtain,  for  an  appropriate  constant  c" 

3/4-8. 
88~ 


i__3/4-8_ 

T{m,n)  <  [  Dc"   '"^  gg     '  +  d  \m'^'^-''n'^'^^^^  log2« 
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Thus,  if  we  choose  r  sufficiently  large  so  as  to  make  the  expression  in  the  brackets 
less  than  ":r+'^.  and  choose  D=2d,  we  obtain 

At 

thus  establishing  the  asserted  inequality,  n 

Remark:  The  preceding  theorem  implies  that  T(m,n)  =  0(m^"*~^n^"*'^^^)  for  any 
8  >  0,  provided  neither  m  nor  n  is  too  small.  This  is  significant  because,  as  shown 
in  [EH]  and  noted  in  section  2,  the  worst  case  total  complexity  of  the  m  cells  con- 
taining the  given  points  can  be  Cl(m^''^n),  which  is  greater  than  0{m^'^n^'^)  for  any 
m  =  o(n^).  In  other  words,  we  have  shown  that  finding  just  the  plane  lying 
immediately  below  each  of  m  given  points  in  an  arrangement  of  n  planes  can  be 
accomplished  faster  than  calculating  the  entire  cell  around  each  point.  While  this 
result  seems  intuitively  plausible,  it  is  noteworthy  that  in  the  2-d  case  no  algorithm 
is  known  as  yet  for  this  seemingly  simpler  task,  which  is  asymptotically  faster  than 
the  algorithm  for  calculating  the  entire  cells  (see  [EGS],  [EGHSSSW]). 

We  next  consider  the  other  problems  Hsted  at  the  beginning  of  the  section.  The 
solutions  to  these  problems  follow  closely  that  of  problem  (a)  given  above,  and 
differ  from  it  only  in  three  key  substeps:  processing  the  data  at  the  bottom  of  recur- 
sion, processing  the  sets  Bv,  and  combining  output  from  children  of  an  internal  node 
V  to  obtain  the  desired  output  for  v. 

Consider  problem  (b),  in  which  we  want  to  count  the  number  of  incidences 
between  m  given  points  and  n  given  planes.  For  a  node  v  at  the  bottom  of  the  recur- 
sion, we  calculate  the  entire  arrangement  of  the  corresponding  n^  planes,  assign  to 
each  (0,1  or  2-dimensional)  face  of  the  arrangement  the  number  of  planes  contain- 
ing that  face,  and  then  locate  in  the  arrangement  the  /My  given  points.  If  a  point  hap- 
pens to  lie  on  a  low-dimensional  face  /,  we  add  to  the  total  incidence  count  the 
number  of  planes  containing  /.  As  above,  all  this  can  be  accomplished  in  time 
0{my  log^/iv). 

When  processing  a  set  fiy  associated  with  some  node  v,  we  calculate  the  upper 
and  lower  envelopes  of  the  riy  planes  in  Ilv,  and  project  them  onto  the  ry-plane  to 
obtain  two  planar  maps  A/"*",  M~ .  We  also  maintain,  for  each  vertex  or  edge  of 
these  maps,  the  number  of  planes  in  Ilv  passing  through  that  vertex  or  edge.  We 
then  project  each  point  p  in  B^,  onto  the  ry-plane  and  locate  its  projection  p*  in,  say, 
M  .  We  can  then  easily  determine  whether  p  lies  on  the  corresponding  upper 
envelope  U  or  not.  If  not,  p  (coupled  with  the  planes  in  11^)  contributes  nothing  to 
the  total  incidence  count.  Otherwise,  if  p  lies  on  a  2-d  face  of  U,  it  contributes  1  to 
the  total  count;  if  it  lies  on  an  edge  or  coincides  with  a  vertex  of  U,  we  add  to  the 
total  count  the  number  of  planes  in  Ily  passing  through  that  edge  or  vertex.  We 
repeat  this  procedure  for  M~ .  It  is  clear  that  this  processing  can  be  accomplished  in 
time  0{(mv  +  nv)  log  riy). 

Finally,  combining  the  output  from  children  of  an  internal  node  v  is  trivial.  We 
simply  add  up  the  incidence  counts  obtained  recursively  for  each  child  iv  of  v,  as 
well  as  the  incidence  counts  for  points  in  the  sets  B^  associated  with  these  children. 
This  takes  only  constant  time. 
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It  can  now  be  easily  checked  that  the  time  complexity  of  the  resulting  algorithm 
obeys  the  same  recurrence  as  the  algorithm  for  problem  (a).  Thus  problem  (b)  can 
be  solved  within  the  same  (randomized)  time  complexity  given  in  Theorem  4.1. 
That  is, 

Theorem  4.2:  One  can  count  the  total  number  of  incidences  between  m  given  points 
and  n  given  planes  in  three  dimensional  space,  in  randomized  expected  time 
Oi[m^''^~^n^"^'^^^  +  m  ]  log^/j  +  n  log  n  log  m),  for  any  8>0. 

Problem  (c)  is  handled  in  a  way  similar  to  that  of  problem  (b),  except  that  now 
we  are  not  concerned  with  the  time  complexity  of  the  algorithm,  but  rather  with 
estimating  the  incidence  count.  In  addition,  the  assumption  that  no  three  points  are 
collinear  turns  out  to  be  quite  strong,  and  yields  sharper  bounds  than  those  obtained 
for  the  previous  problems. 

As  above,  we  construct  a  customized  partition  tree  for  the  dual  points  and 
planes,  with  the  difference  that  the  construction  of  the  tree  stops  this  time  at  nodes  v 
for  which  niy^nl.  At  each  such  node  v  consider  the  arrangement  A  =  A(Uy)  of  the 
«v  planes  in  Ily.  If  a  point  p  that  reaches  v  does  not  lie  on  any  of  these  planes  it  con- 
tributes nothing  to  the  total  incidence  count.  If  it  lies  in  (the  relative  interior  of)  a 
2-d  face  of  A  then  it  contributes  just  1  to  that  count.  Finally,  there  are  0{nl)  inter- 
section lines  containing  all  vertices  and  edges  of  A,  and  none  of  them  contains  more 
than  two  of  the  given  points.  If  such  a  line  /  lies  in  k  of  the  planes  of  IIv,  it  contri- 
butes at  most  2k  to  the  incidence  count,  but  then  we  can  regard  /as    L)  distinct 


lines,  each  formed  by  the  intersection  of  two  of  these  planes.  This  (rather  standard) 
argument  clearly  implies  that  the  total  incidence  count  at  v  is  at  most 
Oinl+niv)  =  Ointy). 

In  processing  the  set  Sv  of  points  whose  dual  planes  have  missed  the  set  11?  for 
some  node  v,  we  consider,  as  above,  the  upper  and  lower  envelopes  of  the  riy  planes 
in  IIv,  and  project  them  onto  the  xy-plane  to  obtain  two  planar  maps  M'*' ,  M~ . 
Each  point  p  in  B^  lies  either  on  or  above  the  upper  envelope  A/"*" ,  or  on  or  below 
the  lower  envelope  M~ .  Clearly  it  suffices  to  consider  only  points  lying  on  these 
envelopes.  If  such  a  point  p  lies  in  (the  relative  interior  of)  some  face  of,  say  M'^ ,  it 
contributes  just  1  to  the  total  incidence  count.  On  the  other  hand,  these  envelopes 
have  0{nv)  edges  (even  when  each  edge  is  counted  with  multiplicity  equal  to  the 
number  of  planes  containing  it),  and  each  can  contain  at  most  two  points.  Thus  the 
overall  contribution  of  By  to  the  incidence  count  is  0(Z7v  +  iv).  where  b^  is  the  cardi- 
nality of  By. 

Finally,  obtaining  the  incidence  count  at  a  node  v  from  the  counts  at  children  of 
V  is  trivially  done  by  summing  up  these  incidence  counts,  plus  the  contributions  to 
this  count  by  points  in  the  sets  B^,  for  all  children  w  of  v.  Hence  if  I{m,n)  denotes 
the  maximum  number  of  incidences  between  m  points  and  n  planes,  with  no  three 
points  collinear,  we  obtain  the  following  recurrence  for  [(m,n). 


0(m)  if  m  >  n^ 

M 

^  I{mi,ni)  +  Oim+n)  if  m  <  n^ 

1  =  1 


lim.n)  = 
where  we  can  assume  that  the  m,,  n,  and  M  satisfy  conditions  (a)-(c)  given  above. 
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By  modifying  the  proof  of  Theorem  4.1,  one  easily  obtains 

Theorem  4.3:  The  maximum  number  of  incidences  between  m  points  and  n  planes, 
provided  no  three  points  are  collinear,  is  0{m^'^~^n'^'^'*'^^  +  m  +  nlogm),  for 
any  8>0. 

It  is  interesting  to  relate  the  above  bound  to  the  0(m^'^n'^'^  +  m  +  n)  bound 
for  the  number  of  incidences  between  in  points  and  n  circles  in  the  plane  [CEGSW]. 
Mapping  points  in  2-d  into  points  in  3-d  and  circles  in  2-d  to  planes  in  3-d  gives  an 
incidence  problem  for  m  points  (no  three  collinear)  and  n  planes,  which  shows  that 
the  bound  in  Theorem  4.3  also  applies  to  incidences  between  points  and  circles  in 
2-d.  As  can  be  seen,  this  bound  is,  however,  slightly  weaker  than  the  bound  of 
[CEGSW]  but  then  it  is  for  a  more  general  problem  (note  that  in  the  above  transfor- 
mation to  3-d,  all  resulting  points  lie  on  the  paraboloid  z  =x^+y^). 

Remark:  If,  instead  of  assuming  that  no  three  points  are  collinear,  we  require  that 
no  three  planes  pass  through  a  common  line,  we  can  use  duality  to  transform  the 
problem  to  that  just  studied,  resulting  in  the  following 

Corollary  4.4:  The  maximum  number  of  incidences  between  m  points  and  n  planes, 
provided  no  three  planes  pass  through  a  common  line,  is 
^(^4/5+28^3/5-S  +  ^  log  „  +  „),  for  any  8>0. 

5.   The  complexity  of  many  cells  in  an  arrangement  of  hyperplanes 

In  this  section  we  extend  the  results  of  section  2  to  arrangements  of  hyper- 
planes in  d  dimensions,  for  any  d>3,  obtaining  upper  bounds  on  the  complexity  of 
m  distinct  cells  in  arrangements  of  n  hyperplanes.  Here  the  complexity  of  a  cell  is 
defined  as  the  number  of  facets  (i.e.  ((i  — l)-dimensional  faces)  bounding  the  cell. 
The  bounds  that  we  obtain  again  almost  match  the  lower  bounds  given  in  [E]. 
Specifically,  we  show 

Theorem  5.1:  The  complexity  of  m  distinct  cells  in  an  arrangement  of  n  hyperplanes 
in  d  dimensions  is  0(m^'^n'^'^  log  n  +  n''~^). 

Proof:  The  proof  proceeds  by  induction  on  the  dimension  d,  but  is  otherwise  quite 
similar  to  the  analysis  in  section  2.  The  base  case  d  =  3  has  been  established  in  sec- 
tion 2.  Suppose  that  d>3  and  that  the  claim  has  been  estabHshed  for  all  d'<d.   Let 

n  =  {iTi Ttn)   be    a    collection    of  n   hyperplanes    in    d  dimensions,    and    let 

Pi Pm  be  m  points  lying  in  m  distinct  cells  of  the  arrangement  A=AiU)  of 

these  hyperplanes.  As  before,  our  task  is  to  bound  the  complexity  of  the  cells  of 
A  (11)  that  contain  these  points. 

Again,  we  partition  11  into  two  subsets  Hi,  112,  each  containing  roughly  n/2 
hyperplanes,  and  refer  to  the  hyperplanes  in  these  subsets  as  the  red  hyperplanes 
and  the  blue  hyperplanes  respectively.  We  obtain  recursively  the  cells  containing 
the  points  p,-  in  each  of  the  subarrangements  A(Ui),  A (112),  and  then  intersect  these 
cells  to  obtain  the  desired  cells  of  A  (11).  Again,  if  a  point  p  lies  in  a  red  cell  R  and 
in  a  blue  cell  B,  then  the  final  cell  E  containing  p  is  simply  R  H  B,  and,  since  all 
these  cells  are  convex  polyhedra,  the  complexity  of  E  is  bounded  by  the  sum  of  the 
complexities  of  R  and  of  B  (in  d>3  dimensions,  this  argument  holds  only  for  the 
number  of  (J  — l)-dimensional  faces  which  is  our  complexity  measure).  Since  R  (or 
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B)  might  contain  more  than  one  point  p,-,  we  again  need  to  establish  a  combination 
lemma  that  controls  the  complexity  of  the  "merged"  cells. 

5.1.   The  Combination  Lemma  for  Hyperplanes 

The  lemma  in  question  is 

Lemma  5.2  (Combination  Lemma  for  Arrangements  of  Hyperplanes):  Given  an 
arrangement  of  /i^  red  hyperplanes  in  J-dimensional  space,  another  arrangement  of 
rib  blue  hyperplanes,  and  m  points  not  lying  on  any  of  these  hyperplanes,  the  total 
number  of  facets  bounding  the  "purple"  cells  containing  these  points  in  the  arrange- 
ment formed  by  all  n  =  rtr  +  rib  hyperplanes,  is  at  most 

p  +  p  +  Oim^'^n'^'^  log  n  +  n''-'^  +  m) 

where  p  (resp.  3)  is  the  total  number  of  facets  bounding  the  cells  of  the  red  (resp. 
blue)  arrangement  containing  these  points. 

Proof:  Again  we  proceed  by  induction  on  d,  with  the  base  case  d  =  3  established  in 
section  2.  Let  pi,  .  .  .  ,pm  be  m  points  in  ^/-dimensional  space.  We  are  given  two 
arrangements  of  hyperplanes,  the  red  arrangement,  consisting  of  tir  red  hyper- 
planes, and  the  blue  arrangement,  consisting  of  rib  blue  hyperplanes.  Let 
^1,  ...  ,^,  be  the  cells  in  the  red  arrangement  containing  the  points  pi  (where  a 
cell  can  contain  more  than  one  point),  and  let  5i,  ...  ,5,  be  the  blue  cells  contain- 
ing these  points.  For  each  p,  let  E,  =  Rs^  n  B^  be  the  "purple"  cell  containing  p,- 
(where  /?,,,  5,,  are  the  red  and  blue  cells  containing  the  point  p,).  By  considering 
only  a  subset  of  the  points  p,-,  if  necessary,  we  can  assume,  without  loss  of  general- 
ity, that  all  the  cells  E,-  are  distinct  (and  thus  have  pairwise  disjoint  interiors). 

As  noted  above,  it  is  easily  verified  that  if  each  red  cell  and  each  blue  cell  con- 
tains just  a  single  point,  then  the  desired  purple  complexity  is  at  most  p  +  p.  Thus,  as 
in  the  2-D  and  3-D  cases,  the  main  step  in  bounding  the  purple  complexity  is  the 
partitioning  of  each  red  cell  that  contains  more  than  one  point  into  subcells,  each 
containing  only  a  single  point,  using  (portions  of)  the  blue  hyperplanes,  and  the 
symmetric  partitioning  of  blue  cells.  This  is  done  using  "red-blue"  cells  and  "blue- 
red"  cells,  in  analogy  to  the  3-D  case.  Our  task  is  to  estimate  the  "red  complexity" 
of  the  red-blue  cells,  namely  the  number  of  red  (sub)facets  bounding  these  cells  (so 
the  blue  facets  along  their  boundaries  are  ignored);  more  precisely,  we  want  to  esti- 
mate the  number  of  additional  red  subfacets  generated  by  the  intersections  of  each 
blue  hyperplane  with  the  current  set  of  red-blue  cells.  Symmetrically,  we  want  to 
estimate  the  additional  blue  complexity  of  the  blue-red  cells. 

Consider  the  cell  splitting  caused  by  the  /-th  blue  hyperplane  tt,-.  Let  Rj  be  a 
red  cell  that  has  already  been  split  into  a  number  of  red-blue  subcells,  and  suppose 
that  IT,  splits  k  =  kj  of  these  subcells  (in  the  sense  of  splitting  the  sets  of  points  they 
contain).  If  ir,  intersects  a  red  or  red-blue  cell  but  does  not  split  the  set  of  points  in 
this  cell  we  still  draw  the  intersection  of  -ir,  with  this  cell,  but  we  can  ignore  its  influ- 
ence on  the  red  complexity  of  this  cell,  since  this  complexity  can  only  decrease  if  it 
changes  at  all.  Since  the  red-blue  subcells  are  convex  polyhedra  by  construction,  the 
intersection  of  ir,  with  each  of  them  is  a  convex  (^- l)-polyhedron  q  (where  these 
polyhedra  have  pairwise  disjoint  interiors),  contained  in  the  intersection 
c  =  Cj  =  Tii  n  Rj  which  itself  is  a  convex  (^- l)-polyhedron.  To  each  facet /of  a  red- 
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blue  subcell  which  is  split  into  two  subfacets  by  it,-,  there  corresponds  a  (d—2)- 
subfacet  of  some  facet  of  c  (which  is  the  intersection  of /with  it,).  Thus  it  follows 
easily  that  the  total  increase  /  in  the  red  complexity  of  the  red-blue  subcells  within 
Rj  as  caused  by  tt,-  is  equal  to  the  number  of  "red"  facets  of  the  red-blue  intersec- 
tion polyhedra  q,  namely  facets  that  lie  along  the  boundary  of  c. 

Let  Af  be  the  arrangement  of  the  rir  (J  — 2)-flats  formed  by  intersecting  it,-  with 
each  of  the  rir  red  hyperplanes,  and  let  At  be  the  arrangement  of  the  /  —  I  (d  —  2)- 
flats  formed  by  intersecting  ir,-  with  the  preceding  /  —  I  blue  hyperplanes.  Let  ^j 
denote  the  total  number  of  red-blue  subcells,  in  all  red  cells,  that  have  been  split  by 
■Hi,  and  let  /,•  denote  the  increase  in  the  red  complexity  of  these  red-blue  cells,  as 
effected  by  it,.  Then  the  preceding  argument  implies  that  /,•  is  equal  to  the  "red 
complexity"  of  |,-  cells  in  the  union  (i.e.  overlay)  of  the  two  arrangements  Af  and 
At ,  where  these  cells  are  the  intersections  of  tt,-  with  the  red-blue  cells  it  splits.  But 
the  red  complexity  of  these  cells  is  bounded  by  their  total  complexity.  Since  this 
arrangement  has  at  most  n  (c?  — 2)-flats,  it  follows,  by  induction  hypothesis  on 
Theorem  5.1  in  d  —  l  dimensions,  that  the  complexity  in  question  is  at  most 

Summing  these  quantities  over  all  rib  blue  hyperplanes,  the  total  increase  in  the  red 
complexity  is  therefore 


O 

which,  as  in  section  2,  is 


"b  "b 

(2  e')-n^'-'^''  log  n+n''-'+  2  ^i 
.=1  1=1     , 


0(;„2/3„d/3  iQg^    ^   ^d-\    ^  ^) 


We  now  obtain  in  a  completely  symmetric  manner  a  partitioning  of  the  blue 
cells  into  a  collection  of  "blue-red"  subcells,  each  containing  just  a  single  point,  and 
derive  a  symmetric  bound  on  their  overall  "blue"  complexity.  Arguing  as  in  section 
2,  we  obtain  the  assertion  of  the  combination  lemma,  n 


5.2.   The  Complexity  of  Many  Cells 

We  now  continue  the  proof  of  Theorem  5.1.  Let  K^'^^{m,n)  denote  the  max- 
imum complexity  of  m  distinct  cells  in  an  arrangement  of  n  hyperplanes  in  d  dimen- 
sions. Applying  the  Combination  Lemma  to  the  two  subarrangements  A{Yii), 
A  (112),  we  readily  obtain  a  recurrence  of  the  form 

K^''Hm,n)  =  2K^'^^{m,nl2)  +  0{m'^'^n'^'^  log  «  +  n'^-'^  +  m)  . 

Again,  we  do  not  continue  this  recursive  process  all  the  way  down.  When  m 
becomes  greater  than  n'^  we  use  the  trivial  bound  K^'^Hm.n)  =  Oin'^)  =  0(m). 
When  this  happens  we  have  m^(2n)'';  it  follows  that 
m  =  m^'^m^'^  =  Oim'^'^n'^'^).  By  the  same  argument,  we  can  drop  the  term  0(m) 
in  the  above  recurrence,  since  we  apply  it  only  for  m^n'^.  It  is  now  easy  to  verify 
that,  since  d/2  >  1,  the  solution  to  the  above  recurrence  is 


K^'^Hm,n)  =  Oim^'^'"^  log  n  +  n'^-^) 
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This  completes  the  inductive  proof  of  Theorem  5.1.  □ 

Corollary  5.3:  The  maximum  number  of  incidences  between  m  points  and  n  hyper- 
planes  in  d  dimensions,  where  it  is  assumed  that  each  point  is  a  vertex  of  the 
arrangement  of  these  hyperplanes,  is  O (m'^'^ n'^'^  log  n  +  n'^"^). 

Proof:  Identical  to  that  of  Corollary  2.3.  a 

Remarks: 

(1)  The  bound  given  in  Corollary  5.3  almost  matches  the  lower  bound 
nim'^'^n'^'^  +  n''"^)  given  in  [E].  It  would  be  nice  to  get  rid  of  the  log  n  factor  in 
our  upper  bound;  we  note  that  it  suffices  to  do  so  in  three  dimensions,  since  in 
higher  dimensions  the  recurrence  for  K^''Hm,n)  will  not  introduce  this  factor. 

(2)  An  obvious  open  problem  is  whether  the  bound  in  Theorem  5.1  also  applies  to 
the  number  of  faces  of  any  dimension  bounding  the  given  m  cells.  Recent  results  in 
[E2]  lend  hope  to  an  affirmative  solution  to  this  problem. 

(3)  Again,  can  this  analysis  be  modified  to  yield  an  algorithm  which  actually  calcu- 
lates the  facets  bounding  m  cells  in  an  arrangement  of  n  hyperplanes  in  d>3  dimen- 
sions? This  seems  to  be  a  considerably  more  difficult  problem  than  the  correspond- 
ing problem  for  three  dimensions,  because  it  calls,  among  other  things,  for  a  pro- 
cedure which  produces  the  facets  bounding  the  intersection  of  two  convex  polyhedra 
in  time  proportional  to  the  number  of  their  facets,  and  no  such  procedure  seems  to 
be  available. 

(4)  Note  that  our  divide  and  conquer  approach  can  be  made  into  an  efficient  algo- 
rithm for  calculating  all  incidences  between  m  points  and  n  hyperplanes.  The  general 
recursive  step  is  trivial,  since  all  we  need  to  do  is  to  sum  up  the  incidences  calcu- 
lated in  the  two  subproblems.  At  the  bottom  of  the  recursion,  when  m^n'^ ,  we  need 
a  procedure  that  finds  the  incidences  between  m  points  and  n  hyperplanes  in  time 
near  0(n'^).  This  can  be  accomplished  using  Clarkson's  randomized  point  location 
algorithm  in  arrangements  of  hyperplanes  (see  [C12]). 
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